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Abstract. We discuss the idea of a "family of L-functions" and describe various methods 
which have been used to make predictions about L- function famihes. The methods involve 
a mixture of random matrix theory and heuristics from number theory. Particular attention 
is paid to families of elliptic curve L- functions. We describe two random matrix models for 
elliptic curve families: the Independent Model and the Interaction Model. 



1. Introduction 

Using ensembles of random matrices to model the statistical properties of a family of L- 
functions has led to a wealth of interesting conjectures and results in number theory. In this 
paper we survey recent results in the hopes of conveying our best current answers to these 
questions: 

(1) What is a family of L-f unctions? 

(2) How do we model a family of L-functions? 

(3) What properties of the family can the model predict? 

In the remainder of this section we briefly review some commonly studied families and 
describe some of the properties which have been modeled using ideas from random matrix 
theory. In Section |21 we provide a definition of "family of L-functions" which has been 
successful in permitting precise conjectures, and we briefly describe how to model such a 
family. In Section El we discuss families of elliptic curve L-functions and show that there 
is an additional subtlety which requires us to slightly broaden the class of random matrix 
models we use. Then in Section E] we discuss how to go beyond the leading-order terms 
which random matrix theory can model, and how one can avoid using random matrix theory 
when modeling a family of L-functions. 

If you only care about elliptic curves and their L-functions you can safely skip to Section |21 

I thank Brian Conrey, Nina Snaith, Matt Young, and Steven J. Miller for many helpful 
conversations. 



1.1. A quick history of families. The idea that collectively the zeros of a single L-function 
behave in a manner that can be modeled statistically started with Montgomery's jMonj work 
on the pair correlation of zeros of the zeta-function. Combined with the large-scale numerical 
calculations of Odlyzko |Odlj . this provided convincing evidence that, to leading order, the 
local statistics of the zeros of the zeta-function, suitably rescaled, were the same as those of 
large random unitary matrices. 

Work supported by the American Institute of Mathematics and by the Focused Research Group grant 
(0244660) from the NSF. This paper is an expanded version of a talk given at the workshop "Random Matrix 
Theory and Elliptic Curves" held at the Isaac Newton Institute, February 2004. 
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A similar collective behavior was noted long ago and is termed the "g-analogue" for Dirich- 
let L-functions. That is, results about all Dirichlet L-functions mod q look just like results 
for the Riemann C-function in t- aspect. For example, the formulas for moments of |C(| + I 
and |L(s,x)| are identical but for replacing t by q. Another early example of collective 
behavior is the pair correlation of zeros of quadratic Dirichlet L-functions |()Slj . Clearly 
something interesting is going on. 

The idea of a family of L-functions with an associated symmetry type began with the 
work of Katz and Sarnak |KSaj . They consider families of function field L-functions, where 
in this case a "family" of L-functions is the set of L-functions associated to a set of curves 
having certain properties. Here the collection of curves must be "natural" in the sense that 
the monodromy group of the family ties it all together. They show that to leading order the 
statistics of the (normalized) zeros of these L-functions, when averaged over the family, are 
the same as the statistics of the (normalized) eigenvalues of random matrices chosen from 
a classical compact group. Here the matrices are chosen uniformly with respect to Haar 
measure, and the size of the matrices scales with the conductor of the L-function. 

For global L-functions there does not (yet?) exist an analogue of "monodromy", but it 
still has been found that to naturally occurring families one can associate a classical group 
of matrices. The zeros of the L-functions have, to leading order, the same statistics as the 
eigenvalues of a randomly chosen matrix from the group. And the appropriately rescaled 
critical values of the L-functions have, to leading order, the same distribution as the "critical 
values" of the characteristic polynomials of the matrices from the group, chosen uniformly 
with respect to Haar measure. The correspondence between the matrices and the L-functions 
involves equating the eigenvalue spacing with the zero spacing, or equivalently setting the 
matrix size equal to the conductor of the L-function. See Section IT^ for more details. 

After the work of Katz and Sarnak there quickly appeared many examples of L-functions 
families behaving in a manner predicted by random matrix theory. Some of the families 
considered were: L-functions associated to holomorphic cusp forms (in either weight or level 
aspect); Dirichlet L-functions (either all or quadratic); and various twists or symmetric 
powers of L-functions. Low-lying zeros were considered by Iwaniec, Luo, and Sarnak |ILSj . 
Rubinstein fRub] , Ozliik and Snyder |US2j . and others. Moments were considered by Iwaniec 
and Sarnak [ISJ, Kowalski, Michel, and VanderKam |iKMVj . Soundararajan S , and others. 

There is an important distinction between the predictions for zeros of L-functions and 
the predictions for moments. For the zeros there is a natural way to normalize: rescale so 
that the average spacing is 1. This rescaling involves the conductor of the L-function and 
the degree of the characteristic polynomial, and this is the source of the principle that one 
chooses the size of the random matrix to equal the (logarithm of the) conductor of the L- 
function. With these normalizations one obtains accurate predictions for the leading-order 
behavior of statistics of zeros of L-functions. However, in addition to its zeros a polynomial is 
determined by an overall scale factor. So one may hope to use the characteristic polynomials 
to model the L-functions, but there will be a correction factor that does not come from 
random matrix theory. The use of the characteristic polynomial to model the L-function 
was begun by Keating and Snaith |KSlj . The first situation which was well understood is 
the moments of L-functions, for which there are explicit predictions for the arithmetic scale 
factor. See p^lTTFlIT^ . 
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2. What is a family? 

There does not yet exist an adequate definition of "family of L-functions" . An attempt 
is made in |(yFKRS] to define a family axiomatically, and we will describe that definition 
here. In that definition the axioms are chosen so that it is possible to produce a plausible 
conjecture for the critical moments of the family. 

Complete details and many examples are in |("FKRSj . so we will just highlight the key 
features of that definition of a family. The main idea is that one starts with a fixed L-function 
and a family of "characters" , and the family of L-functions is produced by twisting the fixed 
L-function by the family of characters. Note that here the term character is used to cover 
more general classes of functions than just Dirichlet characters. 

2.1. L-functions. We wish to define a "family of L-functions", so first we have to give the 
definition we will use for "L-function" . The definition of an L-function which we give below 
is slightly different than what is known as the "Selberg class," but it is conjectured that the 
two are in fact equal. 

Let s = a + it with a and t real. An L-function is a Dirichlet series 

oo 



n=l 



with a„ = Oe{n^) for every e > 0, which has three additional properties. 

Analytic continuation: L{s) continues to a meromorphic function of finite order with at most 
finitely many poles, and all poles are located on the a = 1 line. 

Functional equation: There is a number e with |e| = 1, and a function 7l(s) of the form 

w 

(2.2) ^^^s) = P{s)Q^\[V{\s + ^^,), 

i=i 

where Q > 0, > 0, and P is a polynomial whose only zeros in a > are at the poles of 
L(s), such that 

(2.3) Us) ■■= 1l{s)L{s) 
is entire, and 

(2.4) Us)=eTL{l-s), 



where ii^s) := ^l(s) and s denotes the complex conjugate of s . 

The number w is called the degree of the L-function. That number will also appear in the 
Euler product. 

Euler product: For cr > 1 we have 



(2.5) L{s) = \[\[{l-^,,p 

p i=i 

where the product is over the primes p, and each I'jpjl equals 1 or 0. 

Note that L(s) = 1 is the only constant L-function, the set of L-functions is closed under 
products, and if L(s) is an L-function then so is L{s + iy) for any real y. An L-function is 
called primitive if it cannot be written as a nontrivial product of L-functions. Throughout 
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this paper we assume all L-functions are primitive, although we usually omit the word 
"primitive." 

Conductor: Associated to an L-function is its conductor, a number which measures the "size" 
of the L-function. The paper |CFKRS| introduced a refined notion of conductor which, to 
leading order, is the logarithm of the usual notion of conductor. The refined conductor is 
necessary in order to have any hope of conjecturing the full main term in a general mean 
value of the L-function. Write the functional equation in asymmetric form: 

(2.6) L{s)=eXL{s)L{l-s), 

where Xl{s) = -^—^-r-r-^- Then the refined conductor of L(s), denoted c(L), is given by 
c{L) = \XX^)\. 

2.2. Families of characters. By a family of characters we mean a collection of arithmetic 
functions JF, where each / G is a sequence /(I) = 1, /(2) = 02,/, /(3) = a^j, . . . whose 
generating function 

00 V 

(2.7) ^/(^)=E^=nn(i-/?-^"r^ 

n=l p j=l 

is a (primitive) L-function such that the collection {Lj : / G J^} has some nice properties. 
If we order the L-functions Lj by conductor c(/), then the data {Q; ^i, . . . , n^} in the 
functional equation of Lj should be monotonic functions of the conductor, and the counting 
function M{X) := ^{f G JF | c(/) < X} should be nice. The final condition on the family 
of characters is the existence of an orthogonality relation among the / G JF. Specifically, we 
require that if mi , . . . , mk are integers then the average 

(2.8) 6eim,,...,mk):= lim M{X)-' V /(rm) . . . /(m,)/(m,+i) . . . /(m^) 

exist and be multiplicative. That is, if (mim2 . . . mk, ?^-l?^2 • • • "^fc) = 1, then 

(2.9) 6i{mini,m2n2, . . ■,mknk) = Si{mi, . . .,mk)Si{ni, . . .,nk). 
See Section 3.1 of |CFKRSj for more details. 

2.3. Families of L-functions. Now we create a family of L-functions by starting with a 
fixed L-function 

00 w 



(2.10) M-)=E^=nn(i-7..^"r^ 

n=l p j=l 

Then the elements of our L-function family /^(JF) are the Rankin- Selberg convolutions 



c{s, /) = L/xg(s) = n n n^^ ~ /^p^^tpj-p" 

p i=i j=i 



00 

a 



(2.11) =E 



njxg 



n=l 
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(There may be some issues with the local factors at the bad primes). Note that if it; = 1 or 
V = 1 then 

oo 

(2.12) A^^/) = E^^- 

n=l 

And in particular if Lg is the Riemann zeta-function, then C{s, f) = Lf{s). 

The point of this definition of "family" is that the axioms provide the necessary ingredients 
to apply the recipe in jCJFKRSj to conjecture the full main term in the shifted Kth moment 



(2.13) M{X)-' J2 n LC2 + c^k), 

c{f)<X l<k<K 

or more generally a shifted ratio |CFZj 

(2.14) M(X)- Y. n 

c{f)<Xl<k<K 

Having such a mean value is sufficient to conjecture just about anything you would like to 
know about the zeros and the value distribution of the L-function. See |CSj for examples. 

Note that some families are unions of increasingly large pieces having the same conductor. 
Examples are the Dirichlet L-functions and the L-functions associated to holomorphic cusp 
forms (in either weight or level aspect). For those families it is believed that the heuristics for 
moments will produce a reasonable conjecture for the average over a fixed (large) conductor. 

Although this definition of "family" is useful for certain applications, it lacks the con- 
creteness of the function field case. In particular, there does not yet exist an analogue 
of monodromy for such a family, and computing the symmetry type of the family is not 
straightforward. We discuss this in the next section. 



2.4. Modeling a family of L-functions. Given a family of L-functions one can ask ques- 
tions about its value distribution or about the distribution of its zeros. In most cases current 
technology is not sufficient to answer the interesting questions, so the next hope is to find a 
plausible conjecture. Only recently have such conjectures been found, and the new ingredient 
is to use random matrices to model the family of L-functions. 

The idea is to associate a classical compact group, U{N), Sp{2N), 0{N), S0{2N), or 
S0{2N + 1), to the family. The local statistics of the eigenvalues should agree, to leading 
order, with the corresponding local statistics of the zeros of the L-functions. And, to lead- 
ing order and after compensating by an arithmetic constant, the value distribution of the 
characteristic polynomial 

N 

(2.15) A{z) = Aa{z) = det(/ - A*z) = Y[{l- ^e"*^") 

n=l 

near the point z = 1 should agree with the value distribution of the L-functions near the 
critical point. Here A is an x unitary matrix A and A* is the Hermitian conjugate of A, 
so the eigenvalues of A lie on the unit circle and are denoted by e*^". 
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In the above correspondence the size of the matrix is set equal to the conductor of the 
L-function. (Actually, to an integer close to the conductor, but to leading order such discrep- 
ancies do not matter). To see why this is a natural choice, consider the functional equation 
satisfied by the characteristic polynomial: 

(2.16) Ka{z) = (-l)^det(A)A^,(z-i). 

If we identify (—1)^ det(y4) with e and with Xl{s) then we have a perfect correspondence 
between the functional equations of h.A{z) and L{s), the unit circle playing the role of the 
critical line and z = 1 the critical point. Just as for L-functions, we define the conductor as 
-^z^ evaluated at the critical point, so is the conductor. Note that identifying conductors 
is equivalent to equating the average spacing between the zeros. Values near the critical 
point are modeled using the correspondence A(e~^) ^ L(| + z). 

It remains to identify the matrix group which corresponds to the family. From the func- 
tional equation it is almost possible to determine the group: the only ambiguity is to dis- 
tinguish between S0{2N) and Sp{2N). At one time it was thought that this case could 
be easily resolved because S0{2N) families always arise as "half" of a larger family, the 
other half being modeled by S0(2N + 1). On the other hand Sp{2N) families do not have 
such a "partner". A counterexample to that hope is described in |MD2j . But even if that 
approach were viable, it is unsatisfactory because it relies on the fact that the symmetry 
type can be found among a small list of possibilities. Fortunately, there are other methods. 
One possibility is to compute the 1- and 2-level densities of the family. This usually can be 
done rigorously for functions with small support, and this is sufficient to distinguish among 
the classical compact groups. But again we are relying on the fact that the symmetry type 
can be found on a short list. That objection can be overcome if one can conjecture the level 
densities in the full range, but that can be quite difficult in practice. Another possibility 
is to use the recipe in |CFKRSj to conjecture the moments of the family. This unambigu- 
ously identifies the group, and it also can tell you if the family is not modeled by one of 
those groups. Unfortunately, it is not clear that the recipe in can be applied to all 

interesting families of elliptic curve L-functions, such as the family ^3 given in ()3.6|) . 

2.5. Summary of modeling. Just to be pedantic, we note the following answers to the 
questions posed at the beginning of Section [TJ 

(1) A family of L-functions is a set of L-functions, ordered by conductor, which is built 
in a particular way from a family of characters. The counting function of the family 
should be nice, and the data in the functional equation should be monotonic functions 
of the conductor. 

(2) A family is modeled by associating to it a classical compact matrix group. The 
specific compact group can usually be determined by computing the level densities 
of the low-lying zeros of the family, or by conjecturing the moments of the family. 
The size of the matrices scales with the (logarithmic) conductor of the L-functions. 

(3) To leading order the rescaled local zero statistics of the family are the same as the 
rescaled local eigenvalue statistics of the group. To leading order the critical moments 
of the family equal the critical moments of the characteristic polynomials, up to a 
multiplicative arithmetic constant. For the family we average over the L-functions 
of conductor less than X, and then let X 00. For the matrix groups the averages 
are with respect to Haar measure. 
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The modeling described above will produce leading order asymptotics. To make more 
precise predictions requires other methods, which are described in Section HI 

3. Elliptic curve families 

For the remainder of the paper, E is an elliptic curve over Q, with root number we = ±1, 
and L(s, E) is the L-function associated to E normalized so that s = ^ is the critical point. 
If F is any family of elliptic curves we write F = F^ U F~ where or F~, respectively, are 
the curves E & F with we = +^ or we = We write Ea^b for the curve y"^ = + ax + b. 

Most of the information in this section can be found in recent papers by Steven J. Miller 
and Eduaro Duefiez [Ml IMUTl IMD2J . Nina Snaith [SnTllHn2] and Matthew Young jYTHY2l 
IY3j . The author of this paper is just trying to convey the current understanding of the 
relationship between families of elliptic curve L-functions and random matrix theory: he 
makes no claim to any of the ideas presented here. 

3.1. Families with a given rank. The following question is not well posed: 

What is the correct random matrix model for the L-functions of a family of elliptic curves 
having a prescribed rank r ? 

The question is not well posed because there are (at least) two reasonable models, both 
of which seem to be appropriate for certain families of elliptic curves. We will first examine 
the simplest case of rank r = 1. 

Consider the following families of rational elliptic curves: 

(3.1) J-i(X) = {E,,b : \a\ < \b\ < X^} 
and 

(3.2) J-2(X) = {^,,,2 : |a| < xi, b' < X^}. 

Note that the point (0, b) on Ea^b^ almost always has infinite order, so almost all of the curves 
in J-'2 have rank at least 1. 

Let's consider J-"^ and J-'2 . In both families we have L{^,E) = because We = — 1- 
However, for E E we could have said L{^,E) = because rank{E) > 1. The fact 
that the zero at L{l,E) for E e J^2 was constructed, instead of just arising from parity 
considerations, has a profound influence on the behavior of the L-function near the critical 
point. 

To understand the influence of the critical zero, we first consider the distribution of L'(i, i?), 
which assuming standard conjectures is nonzero for almost all curves in our we = —1 fami- 
lies. We have the following conjectures from |Y2j : 



(3.3) Yl L\\,Ef^c,{k){\ogXf^-^+^y' 
while 

(3.4) Yl L\\,Ef^c2m^ogxr>^-^y' 
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As the formulas show, the behavior at the critical point is different for the two families, 
even though both families could be described as "a rank 1 family of elliptic curves." In 
particular, we see that the derivative L'{^,E) tends to be smaller for E G J-'2 . This can be 
explained by the tendency for the low-lying zeros of L{s, E) to be closer to the critical point 
for E G J^2 ■ That is, should have more low-lying zeros, which will cause the L-function 
to stay small near the critical point, and so its derivative will also be small. To make this 
idea precise we consider the 1-level density of the zeros. 

Let < ^E,i < lE,2 < lE,z < ■ ■ ■ denote the imaginary parts of the zeros of L(s, E) in the 
upper half of the critical strip. Note that we have omitted the zero(s) at the critical point. 
The one-level density of the family F{X) is defined to be the function Wi which satisfies 

(3.5) /'^W^iW^^' 

' ^ E&F{X) j 

as X ^ oo, for nice functions (p. That is, Wi measures the density of the zeros of the family. 

The observation about the relative size of L'{^, E) can be restated as: the one-level density 
for the family should be more concentrated near than the one-level density for the 
family jFj". By using random matrix theory and some other ideas we explain below, it is 
possible to produce a precise conjecture for the one-level densities of these families. These 
are given in Figure ITTl The functions are rescaled so that the average spacing between zeros 
is 1. In the next section we explain where those conjectures came from. 




12345 12345 



Figure 3.1. Conjectured one-level density of the noncritical L-function zeros of 
the family jFf (left) and (I'ight). 



3.2. Two models for two kinds of families. The plots in Figure IHTTl are familiar. The 
plot on the left is the rescaled one-level density of the eigenvalues of matrices from the group 
S0{2N + 1), in the limit as ^ 00. The plot on the right is the rescaled one-level density 
of the eigenvalues of matrices from the group S0{2N), in the limit as N —>■ 00. Given those 
plot, what are our models for the two families? 

Consider the following ways to make a polynomial f{z) which has real coefficients, all its 
zeros on the unit circle, and (almost surely a simple) zero at 2; = 1: 

• The characteristic polynomial of a matrix in S0{2N + 1) 

• {z — 1) times the characteristic polynomial of a matrix in S0{2N) 
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It should be clear that those two examples will have the one-level densities pictured in 
Figure IXTl These examples are the simplest cases of the two most commonly studied higher 
rank families of elliptic curve L-functions, which we now describe. 

Suppose Et is a curve = + a{T)x + b{T) of rank r over Q{T). Consider the following 
two families of rational elliptic curves: 

(3.6) J^siX) = {Ea,b ■■ \a\ < X^, \b\ < X^, rank(^) > r} 
and 

(3.7) ^4(X) = {Et : \a{t)\ < X^ \b{t)\ < X-^ , t E N}. 

As in our rank 1 example, we have J-'a{X) C J-'3{X). Again we consider the subfamilies 
according to the sign of we- Let F be either of the above rank r families. If r is odd then 
almost all the curves in F~ have rank r, and almost all the curves in have rank r + 1. 
If r is even then the F'^ curves have rank r and the F~ curves have rank r + 1. It is 
conjectured that JF^ and JF^ are approximately the same size provided E has at least one 
place of multiplicative reduction. See [HI- If r is even then it is possible that J-'^ and jFg" 
are approximately the same size. 

We describe the two models which are believed to correspond to these families. The names 
for these models was coined by Steven J. Miller. 

3.3. Selecting to have zeros. The Interaction Model 

We are modeling a family that arises by restricting a much larger family to a subfamily 
having at least r zeros at the critical point. A matrix model for this family can be described 
as follows: start with SO{M) where M = 2N or 2N + 1 depending on whether we are 
modeling J-"^ or jFg", and restrict to those matrices having 1 as an eigenvalue of multiplicity 
at least r. That is, you take matrices in S'O(M) and drag r zeros to the critical point. 

There are some bad things about this model. First, it is a set of matrices, but it is not a 
group. And while it is a perfectly well-defined set, it is a measure zero subset of SO{M), so 
there is no canonical way to restrict Haar measure to it. 

One solution, which has been analyzed by Snaith |Snlj and Duenez |D1 IMDlj is to first 
restrict to those matrices which have r eigenvalues in [—£,£], and then let e — > 0. The 
resulting measure is the same as one obtains by taking Haar measure on SO{M), formally 
substituting ^ = for r of the eigenvalues, and then omitting those terms which vanish 
identically. With the eigenvalues given by e*^-'', the induced measure on that set is 

M 

(3.8) C{M,r)l[{l-cose,r J] ( 

cos 9j — cos OkY'dOi ■ ■ ■ dO^ 

j=l l<j<k<M 

where C(M, r) is a normalization constant. 

It is instructive to look at the one-level density for such matrices. The one-level density 
is given by [EnH IMDT] 

(3.9) [jldO-^] + Jldd^) - ^-^Jr-l^O^Ur-liO-^?) ■ 

There is some numerical evidence |MDlj that this model is accurate. For this model it is 
possible to compute the critical moments of the characteristic polynomials |Snlj , but it does 
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Figure 3.2. One-level density for SO{M), restricted to have exactly r eigenval- 
ues at 6* = 0, for r = 1, 2, 3, 4, 5. 

not seem that all the ingredients are available to use the recipe in jCFKRSj to conjecture 
the moments of the family. 

Note that when r = we recover S0{2N) and when r = 1 we have S0{2N + 1). 

3.4. Imposing zeros. The Independent Model 

We are modeling a family that has an rth order zero at the origin which arises from an 
explicit construction. A model for this situation can be found by assuming that the extra 
critical zeros are just inserted at the critical point, and all the other zeros ignore them. That 
is, start with a matrix in SO{M) where M = 2N — r or 2A^-|-1 — r depending on the parity of 
r and the sign of we- Then the polynomial which models the L-function is the characteristic 
polynomial of the matrix, multiplied by {x — ly. The one-level density only depends on the 
parity of r and the sign of we, and will be one of the functions shown in Figure IXTl 

One can phrase the model strictly in terms of matrices by saying that the model is given 
by the group 



where Irxr is the r x r identity matrix and M = 2N — r or 2N -|- 1 — r is chosen according 
to the parity of r and the sign of the functional equation. 

There is numerical evidence |MDlj . level density calculations [Ml IMDll 1511 lYlj and conjec- 
tures for moments (Y2j that this model gives accurate predictions for some specific families. 

3.5. Some issues. It is worth repeating that the above models, even if they are correct, are 
only intended to capture leading-term asymptotics. Computer experiments [MDQ find that 
the low lying zeros of the family JF4 exhibit some anomalous behavior which presumably will 
disappear when larger examples are computed. It is entirely possible that family JF3 is a 
union of families of the form JF4, and this may contribute to a bias in those numerics. In 
families of type JF4 the generators of the set of rational points have very small height, and 
this may also introduce a bias. That is, the heights of the generators are on the order of the 



(3.10) 
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logarithm of the conductor, while it is more typical to have the heights as large as a power 
of the conductor. See Silverman |Si2j . Chapter 10, for a discussion of heights of generators. 

It is possible that more accurate predictions of the one-level density (using methods de- 
scribed in the next section) will show better agreement with the data. Those methods are 
also able to give extremely precise predictions for the moments of the L-functions |Y2j , and 
these give support to the models. 

It is not clear that these families of elliptic curves give rise to families of L-functions as 
described in Section El For the purpose of conjecturing the moments of the family (which 
is why that definition of 'family' was developed), the key property is the orthogonality 
relation ()2.8p . For specific families of type ^4 it should be possible to evaluate such sums. The 
result is likely to be quite complicated, as in |Y2j . A subtle problem is that the parameter X 
in the elliptic curve families is approximately the discriminant, not the conductor. By 
Szpiro's conjecture the logarithm of the discriminant is within a factor of 6 of the logarithm 
of the conductor, so it is possible that ordering by discriminant is almost as good as ordering 
by conductor. It seems reasonable to model by setting iV, the size of the matrix, equal to 
logX, since what else would you choose? If that choice is correct it suggests that X is close 
to the discriminant most of the time. This has been shown for some families |Ylj . 

One can cook up an elliptic curve family which presumably is a hybrid of the models 
described here: take a rank r family of type JF4 and restrict to those curves having rank at 
least r + 2. If one makes the reasonable assumption that the "extra" zeros created by this 
process do not interact with the original r zeros imposed at the critical point, then one can 
use methods similar to CKRS to predict how many curves are in the restricted family. I 
am not advocating a reckless proliferation of elliptic curve models, but merely noting that 
even if the two models described here are correct and can be refined to predict lower order 
terms, they may not cover all families of interest. 

In Figure 13.31 one can see that if r is large then you are unlikely to find noncritical zeros 
close to the critical point. The name "repulsion" has been given to this phenomenon. The 
logic behind the name is that the "lowest zero" is further from the critical point than it would 
be if there were not a multiple critical zero. Unfortunately, the "lowest zero" is not a well 
defined object. If you drag the lowest zero to the critical point then the other zeros follow 
it toward the critical point, and at the moment you increase the order of the critical zero 
there becomes a new "lowest zero" . If there is an rth order critical zero and you count zeros 
correctly, then the "lowest zero" is actually the (r + l)st zero, and it is likely to be closer to 
the critical point than a typical (r + l)st zero. In other words, the word "attraction" more 
accurately describes the situation! There is no reason to change the current terminology, 
but keep in mind that in the model where one restricts to those matrices having multiple 
eigenvalues at 1, the other eigenvalues have actually moved towards the critical point. 

4. Refined modeling 

Random matrix theory is useful for making leading-order asymptotic predictions about 
families of L-function. To understand the finer behavior of the family one must use heuristic 
techniques from number theory. 

There are two main refinements to the leading term asymptotics. First, L-function fam- 
ilies are ordered by conductor, and we use the conductor to determine the appropriate size 
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matrices for our model. There will be a discrepancy between the limiting behavior for large 
matrices and the behavior for finite size matrices. For some quantities, such as the nearest 
neighbor spacing of the zeros/eigenvalues, it is difficult to see the difference between the 
asymptotics of the distribution and the distribution for moderate size matrices. For other 
quantities, such as the value distribution, for any computable range there is a notable differ- 
ence between the limiting quantity and the values that can be computed. See Keating and 
Snaith |KSlj . for a spectacular example concerning the value distribution of 3?logC(| + H)- 
By a theorem of Selberg that is Gaussian in the limit t —>■ oo, but for finite t it differs from a 
Gaussian in the same way as the characteristic polynomial of an appropriately sized random 
unitary matrix. 

The second issue is the fact that there are lower order terms, and the lower order terms for 
L-functions involve arithmetic factors, while the lower order terms of random matrices do not. 
Thus, the general shape of expressions from random matrix theory can reveal what to expect 
for L-functions, but the arithmetic "correction terms" must be determined in some other 
manner. For zero spacings there are no arithmetic corrections in the leading order terms. 
For moments of L-functions the leading order correction terms are fairly straightforward 
to determine |CG| ICF| IKS2j . Just about everything else is quite subtle and one needs 
sophisticated number-theoretic methods in order to make sensible conjectures. 

For moments of L-functions such conjectures are covered in detail in |CFKRSj . Matt 
Young |Y2j used these heuristics to compute the full main term for various families of el- 
liptic curve L-functions. Those conjectures give ()3.3|) and ()3.4|) as special cases. Thus, the 
heuristics appear to correctly handle a variety of interesting families. (For the family JF3 
in ()3.6|) our current understanding of the distribution of the coefficients does not seem 
adequate to conjecture the moments of the family. Even finding the leading order arithmetic 
factor seems difficult in this case. However, random matrix calculations |iSnl|| predict the 
general shape of the moments.) 

For quantities involving zero statistics, moments are insufficient and one needs averages 
of ratios of the L-functions. This is addressed in |CFZj . For example, the expected value of 
the ratio 

(4 1) ^'^ + °-^) 

is sufficient to determine the one-level density of the family L{s,f), including the lower 
order correction terms due to arithmetic effects. See CS for many examples. It should be 
possible to use these methods to conjecture the ratio (j4. ip for various higher rank elliptic 
curve families, and thus give a precise conjecture for the one-level density. As of this writing 
this has not been done, but probably it will have been done by the time this paper appears 
in print. 

In summary, by choosing the matrix group and the size of the matrices appropriately, and 
with the appropriate arithmetic correction factor, a family of L-functions can be modeled 
by the characteristic polynomials of a collection of matrices. The example of elliptic curve 
L-functions where the elliptic curves are selected to have large rank shows that the collection 
of matrices may not be a group. In order to capture the lower order terms one must use 
heuristics from number theory which do not explicitly involve random matrix theory. Those 
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heuristics also recover the leading order behavior which previously required random matrix 
theory. 
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